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A PLASTICITY  FORMULATION 
FOR  CYCLIC  INELASTIC 
STRUCTURAL  ANALYSIS 

1 INTRODUCTION 

Background 

A structure’s  energy  absorption  capacity  is  impor- 
tant to  its  safety  dicing  dynamic  loading  such  as  that 
occurring  during  an  earthquake.  Inelastic  hvsterctic 
deformation  is  a major  source  of  energy  dissipation 
in  structures  under  such  loadings.  Therefore,  esti- 
mates of  energy  dissipation  should  be  based  on  an 
inelastic  deformation  analysis  appropriate  to  cyclic 
loading  conditions.  Such  an  analysis  would  be  useful 
for  evaluating  plastic  structural  design  procedures  as 
well  as  inelastic  instability  and  cumulative  fatigue 
damage  criteria.  Inelastic  analyses  and  design  proce- 
dures will  also  help  achieve  safe  and  economical 
structures  through  effective  use  of  their  strengths  in 
the  inelastic  range. 

A review  of  current  structural  analysis  procedures 
and  plasticity  formulation  indicates  that  the  material 
constitutive  equations  (multiaxial  stress-strain  rela- 
tionships) used  in  these  procedures  do  not  ade- 
quately represent  material  characteristics  under 
cyclic  loading.  Using  such  an  analysis  would  cause 
serious  errors  in  estimating  energy  dissipation, 
ductility,  anil  inelastic  deformation. 

Extensive  low-cycle  fatigue  research  during  the 
past  two  decades  has  provided  considerable  insight 
into  various  plastic  hysteresis  phenomena  and  cyclic 
history  dependence  of  the  uniaxial  stress-strain 
behavior  of  metals.  Based  on  these  findings,  a uni- 
fied approach  of  characterizing  the  cyclic  uniaxial 
properties  has  been  developed.1  Sufficient  experi- 
mental data  on  cyclic  material  characteristics  under 
multi-axial  stress  states  are  not  yet  available.  There- 
fore. the  only  feasible  approach  is  to  develop  a 
formulation  within  the  framework  of  the  established 
classical  theory  of  plasticity,  w hich  can  simulate  the 
salient  cyclic  uniaxial  material  characteristics. 

Purpose 

I he  purpose  of  this  report  is  to  present  a plasticity 
formulation  for  structural  materials  that  adequately 

'Ft.  R.  Jhansalc.  "A  Friction  Stress  Method  tor  the  Cyclic  In- 
el. otic  Behavior  ol  Metals.  I ransactions  of  f hr  3rd  International 
t'onierence  on  Sinn rural  Mechanics  in  Reactor  Technology.  Vol 
5.  Paper  15  4 (1475) 


represents  all  the  important  cyclic-  material  proper- 
ties and  that  ean  be  readily  incorporated  into  exist- 
ing finite  element  nonlinear  structural  analysis  com- 
puter codes. 

Approach 

A previous  CERL  report2  extended  a cyclic  uni- 
axial stress-strain  model  to  the  generalized  state  of 
stress  by  defining  an  equivalence  between  the  two 
states.  One  problem  with  this  approach  was  the 
necessity  of  reducing  a general  nonproportional 
loading  path  into  only  two  events  (loading  and  un- 
loading) to  define  load  reversals. 

This  problem  has  been  eliminated  in  this  report  by 
adopting  a different  approach  based  on  Mroz’s  con- 
cept of  multiple  and  nesting  hardening  surfaces.3 
The  translational  and  dilatational  properties  of  these 
surfaces  are  determined  from  a set  of  cyclic  uniaxial 
properties  developed  by  Jhansale.4  The  inelastic- 
deformation  analysis  using  this  formulation  should 
help  achieve  an  accurate  and  reliable  evaluation  of 
structural  response  under  dynamic  loading. 

2 LITERATURE  REVIEW 

Many  recently  proposed  mathematical  models5'7 
describe  deformation  of  materials  in  the  inelastic- 
range.  However,  for  rate-insensitive  materials,  only 
the  classical  theory  of  plasticity  and  its  various  modi- 
fications have  been  used  successfully  in  general 
structural  analysis.  This  is  due  to  their  relative- 
mathematical  simplicity  and  their  ability  to  simulate- 
actual  material  behavior  at  least  under  monotonic 
loadings.  The  classical  theory  of  plasticity  with  a 
constant  yield  surface  (perfect  plasticity  model)  or  a 
monotonically  expanding  yield  surface  (isotropic 

?J  F.  McNamara  and  S.  K.  Sharma.  Isotropic-Kinematic 
Hardening  Model  for  Elastic-Plastic  Structural  Analysis,  Techni- 
cal Report  M l 48/ A DAO  1 4945  (Construction  Engineering  Re- 
search Laboratory  |CHRL|.  August.  1975). 

JZ.  Mroz,  "An  Attempt  to  Describe  the  Behavior  ol  Metals 
Under  Cyclic  Loads  Using  a More  General  Workhardening 
Model.”  Acta  Mcchanica,  Vol  7,  No.  2-3  (19b9).  pp  199-212. 

4Jhansalc.  "A  Friction  Stress  Method  tor  the  Cyclic  Inelastic 
Behavior  of  Metals." 

SJ.  R Rice.  "On  the  Structure  of  Stress  Strain  Relations  tor 
Time- Dependent  Plastic  Deformation  of  Metals."  Journal  of 
Applied  Mechanics.  Vol  37  (197())_  pp  728-737. 

6T.  Mura.  “Continuum  Theory  of  Dislocations  and  Plasticity.” 
Mechanics  of  Generalized  Continua.  I Kroner,  ed.  (Springer 
Verlag.  I9b8).  pp  2b9  278. 

7K.  C.  Valanis.  "A  Theory  ot  Viscoplasticitv  Without  a Yield 
Surface."  Archives  Mechanics.  Vol  23  (I97|t.  pp  SI  7 SSI 
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hardening  model)  has  been  presented  in  detail.8,9  A 
constant  or  expanding  yield  surface  conflicts  with 
the  experimental  Baushinger  effect  phenomenon, 
which  is  denoted  by  a reduced  yield  strength  after  a 
stress  reversal  following  plastic  straining.  To  simu- 
late the  Baushinger  effect,  I’rager8 * 10  introduced  a 
kinematic  hardening  model  in  which  yield  surface 
translates  without  rotation  or  distortion  in  the  direc- 
tion of  plastic  strain  increments.  This  kinematic 
hardening  model  was  modified  by  Ziegler11  to  make 
translation  of  the  yield  surface  invariant  under  a 
reduction  of  the  dimension  of  stress  space. 

The  use  of  the  classical  plasticity  theory  with 
isotropic  or  kinematic  hardening  models  requires  a 
material  constant  which  is  usually  determined  from 
an  experimental  uniaxial  stress-strain  curve.  A 
correspondence  is  established  between  multiaxial 
and  uniaxial  stress-strain  paths  by  defining  equiva- 
lent increment  of  plastic  strain  from  multiaxial 
plastic  strain  increments.  The  material  constant  is 
then  determined  from  the  slope  of  the  experimental 
uniaxial  stress-strain  curve  at  a point  which  corre- 
sponds to  current  stress  state  for  the  multiaxial  case. 
This  approach  for  determining  the  material  constant 
experiences  several  difficulties  when  cyclic  plastic 
deformations  are  involved.  If  only  a monotonic 
experimental  uniaxial  curve  is  used,  then  multiaxial 
plastic  strains  in  any  direction,  including  the  com- 
pletely reversed  direction,  must  be  reduced  to  a 
monotonically  increasing  plastic  strain.  If  reverse 
plastic  straining  for  the  uniaxial  curve  is  also  per- 
mitted. then  every  increment  of  multiaxial  plastic 
deformation  must  be  defined  as  a case  of  either  load- 
ing or  reversed  loading.  There  is  no  smooth  transi- 
tion from  a case  of  continued  loading  to  a ease  of 
completely  reversed  loading.  Moreover,  an  analytical 
model  for  generating  cyclic  uniaxial  stress-strain 
must  be  formulated.  This  curve  cannot  be  obtained 
experimentally  because  the  history  of  deformation  is 
not  known  a priori. 

A definition  for  determining  loading  or  reversed 
loading  from  multiaxial  plastic  deformation  and  a 
rheological  material  model  for  steel  have  been  used 
in  the  application  of  the  kinematic  hardening 

8W.  Prager  and  P.  G.  Hodge.  Theory  of  Perfectly  Plastic  Solids 
(John  Wiley.  1951). 

’R.  Hill.  The  Mathematical  Theory  of  Plasticity  (Oxford . 1950). 

10W  Prager.  "A  New  Method  of  Analyzing  Stresses  and  Strains 

in  Work  Hardening  Plastic  Solids."  Journal  of  Applied  Meehan 

us.  Vol  21  (1956).  pp  491  496, 

"H.  Ziegler.  "A  Modification  of  Prager' s Hardening  Rule." 

Quarterly  of  Applied  Mathematics.  Vol  17  (1959).  pp  55-65. 


theory.12  A similar  procedure  for  a combined  iso- 
tropic-kinematic hardening  theory  has  been  used 
where  a method  to  determine  another  material  con- 
stant for  isotropic  hardening  is  also  given.13  Re- 
cently, a multiaxial  plastic  deformation  in  any  direc- 
tion has  been  considered  as  a case  of  either  con- 
tinued loading,  partially  reversed  loading,  or  fully 
reversed  loading.14  A cyclic  uniaxial  curve  has  been 
represented  by  branches  of  a Ramberg-Osgood 
curve,15  which  is  essentially  a power  law  relationship 
between  stresses  and  strains.  The  three  material 
parameters  of  this  curve  are  adjusted  according  to  a 
scalar  product  of  previous  and  current  normal 
vectors  in  the  direction  of  deformation. 

A more  consistent  theory  which  can  account  for 
the  directional  dependence  of  reversed  multiaxial 
plastic  deformation  has  been  proposed  by  Iwan16 
and  by  Mroz.17  In  Mroz’s  model,  a uniaxial  stress- 
strain  curve  of  an  initially  isotropic  material  is 
approximated  by  piecewise  linear  segments  of  con- 
stant tangent  plastic  moduli.  In  a stress  space,  this 
approximation  is  represented  by  concentric  and 
similar  hypersurfaces  defining  regions  of  constant 
plastic  tangent  moduli.  The  first  of  these  surfaces  is 
the  conventional  yield  surface.  The  yield  surface 
translates  during  plastic  deformation,  coming  in 
contact  with  and  then  pushing  larger  hypersurfaces. 
Motions  of  these  surfaces  are  governed  by  the  condi- 
tions that  they  preserve  their  orientations  and  that 
no  two  surfaces  can  intersect  each  other.  These  con- 
ditions give  rise  to  a kinematic  hardening  rule  that 
differs  from  those  of  Prager18  and  Ziegler.19 


1!F.  G.  Plummer.  A New  Look  at  Structural  Energy  Dissipa 
lion.  Technical  Report  M-82  AD#780801  (CKRL.  May  1974). 

1!J.  F.  McNamara  and  S.  K Sharma.  fsotropir-fCiriemalic 
Hardening  Model  for  Elastic-Plastic  Structural  Analysis.  Techni- 
cal Report  M 148  A DAO  14945  (CFRL,  August  1975) 

“M.  A.  F.isenbcrg.  "A  Generalization  of  Plastic  Flow  Theory 
With  Application  to  Cyclic  Hardening  and  Softening  Phenome 
na,"  Journal  of  Engineering  Materials  and  Technology.  Vol  98. 
Series  H.  No.  1 (July  1 976).  pp  221  228. 

,8W.  Ramhcrgand  W R Osgood.  Description  of  Stress  Strain 
Curves  hy  Three  Parameters  TN  902  (National  Advisory  Com- 
mittee for  Aeronautics.  1941). 

»W.  D.  Iwan.  "On  a Class  of  Models  for  the  Yielding  Behavior 
of  Continuous  and  Composite  Systems."  Journal  of  Applied 
Mechanics.  Vol  .14  (1967).  pp  612-617. 

17Z.  Mroz.  "On  the  Description  of  Anisotropic  Workharden 
ing."  Journal  of  Mechanics  of  Physical  Solids.  Vol  15  (196"1). 
pp  16.1-175. 

>'W  Prager,  "A  New  Method  of  Analyzing  Stresses  and  Strains 
in  Work-Hardening  Plastic  Solids."  Journal  ot  Applied  Met  han 
ics.  Vol  21  (1956).  pp  491  496 

18H  Ziegler.  "A  Modification  ol  Prager's  Hardening  Rule." 
Quarterly  ot  Applied  Mathematics.  Vol  l"1  11950).  pp  55  65 
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Recently.  Krieg20  and  Datalias  and  Popov21  have 
proposed  a simplification  of  Mroz’s  model  through 
consideration  of  only  two  surfaces — a yield  surface 
w Inch  is  enclosed  by  a limiting  surface.  Both  of  these 
surfaces  arc  allowed  to  translate  and  change  in  si/e 
during  plastic  deformation.  The  rates  of  translation 
and  the  changes  in  si/e  are  assumed  to  be  functions 
of  a vector  from  the  current  stress  state  to  a projected 
stress  state  on  the  limiting  surface.  However, 
Datalias  and  Popov  have  left  these  functions  un- 
defined because  of  a lack  of  experimental  data. 
Kricg  has  defined  these  functions  intuitively  to 
match  experimental  data  for  a uniaxial  specimen, 
but  their  validity  for  multiaxial  deformations  has  not 
been  demonstrated.  Moreover,  the  use  of  only  two 
surfaces  is  inadequate  for  a satisfactory  simulation 
of  the  “memory"  effect22  in  which  a material  appar- 
ently remembers  certain  previous  points  of  stress- 
strain  reversals.  The  material  exhibits  discontinui- 
ties in  its  stress-strain  response  at  these  points  on 
subsequent  loading  or  unloading. 

In  most  of  the  theories  presented  above,  isotropic 
hardening  is  taken  into  account  by  changing  the 
radius  of  the  yield  surface  as  a function  of  a monoto- 
nieallv  increasing  parameter  (e.g..  equivalent  plastic 
strain,  plastic  work).  This  usually  implies  that  the 
radius  of  the  yield  surface  increases  (or  decreases) 
monotonieally,  reaching  a stable  value  asymptoti- 
cally. An  obvious  drawback  of  this  approach  is  that 
the  yield  radius  is  restricted  to  increase  (or  decrease) 
monotonieally  regardless  of  whether  the  amplitude 
of  loading  cycles  is  increased  or  decreased.  Once  the 
yield  radius  has  reached  a stable  value,  it  remains 
unchanged  even  if  the  amplitude  of  loading  cycles  is 
changed  significantly.  To  overcome  this  problem, 
the  idea  of  nesting  hypersurfaces  is  extended  in  the 
formulation  presented  in  this  report  to  define  a field 
of  isotropic  hardening  moduli.  All  material  con- 
stants required  by  the  present  approach  are  obtained 
front  cyclic  uniaxial  experiments  described  briefly  in 
Chapter  .f. 


'”K  I)  Knit  A Practical  two  Surface  Plasticity  Theory." 
Journal  of  \pplt>  J Vtci  hiinii  s Vol  42  1197S),  pp  o4 1 64b 

t Datalias  amt  i I1  Popes.  "A  Model  of  Nonlineartv 
Hardening  Malcrials  lor  C i an  pit’s  I.oading."  Actii  Mrchanicu. 
Vol  21  IIT5I.  pp  1 7,1  |92 

JJII  K Jhaitsaia  amt  I H topper,  An  Fngineering  Analysis 
o!  I he  Inelastic  Stress  Response  ol  a Structural  Metal  Under  Vari- 
able Is  die  Strain.  I ipmmcntatinn  anil  hiithtrr  Pn-iliilinn 
Sll*  5 1 u i American  Sociels  tor  testing  and  Materials.  107!) 


s) 


3 


CYCLIC  UNIAXIAL 
STRESS-STRAIN  BEHAVIOR 


In  order  to  meaningfully  describe  the  cyclic  uni- 
axial material  properties,  a brief  outline  of  the 
plastic  hysteresis  phenomena  of  metals  and  a re- 
cently developed  unified  approach23  for  their  quanti- 
tative characterization  is  necessary. 

Transient  Behavior 

Cycle-dependent  variations  in  the  hystcrctic 
behavior  of  metals  are  usually  observed  under  con- 
stant stress-  or  strain-limits  cycling  conditions. 
These  variations  are  classified  as  cyclic  hardening, 
cyclic  softening,  cyclic  creep,  and  cyclic  relaxation. 
In  general,  these  phenomena  are  transient,  and  a 
saturation  or  stable  state  is  rapidly  approached  at  a 
continuously  decreasing  rate  with  cycles.  In  general, 
materials  exhibit  an  initial  cyclic  hardening  or 
softening,  depending  on  their  initial  condition  and 
sometimes  on  the  magnitude  of  the  stress  or  strain 
amplitude  imposed.  Whenever  the  stress  or  strain 
limits  are  changed,  a transient  behavior  ensues  and  a 
new  saturation  state  corresponding  to  the  new  stress 
or  strain  limits  is  approached. 

Based  on  the  observations  of  a large  number  of 
structural  metals  of  widely  ranging  strengths  and 
cyclic  characteristics,  it  has  been  shown24-25  that 
transient  behavior  is  essentially  caused  by  a history- 
dependent  variation  in  the  material's  yield  strength. 

I his  variation  is  effected  by  substructural  changes 
caused  by  the  reversed  plastic  strains  in  the  material. 

Figure  1 schematically  illustrates  a cyclic  harden- 
ing situation  that  compares  the  first  three  hysteresis 
branches.  Each  hysteresis  branch  is  composed  of  an 
initial  linear  elastic  portion  and  a nonlinear  inelastic 
portion.  In  the  ease  of  cyclic  hardening,  the  initial 
elastic  portion  increases  in  length  from  one  branch 
to  the  next,  but  the  elastic  modulus  remains  tin 
changed.  The  nonlinear  inelastic  portions  remain 
virtually  unchanged  in  shape  during  transient 
behavior.  The  history -dependent  change  in  the 

”H  R Jhansale.  "A  Friction  Stress  Method  lor  the  ( vcltc  In 
clastic  Behavior  ot  Metals."  / run.stu  linns  ol  ihr  Ini  Intrrnaiional 
lonlinin  i on  Structural  Sfci  homes  m Kcactor  lcchnolop\  Vol 
5.  Paper  I S 4 iprsi 

'*11  R Jhansale.  A Friction  Stress  Method  tot  the  ( vein  In 
elastic  Behavior  ol  Metals  " 

n!l  R Jhansale.  "A  New  Parameter  tor  the  Ifvstervtu  Stress 
Strain  Behai  tor  ot  Metals . " Journal  ot  I npmccnnp  \tai*-nah  ttfti/ 
I rt  hnologt  Vol  9’.  No  I IWSI.  pp  It  IK 
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(a)  Three  Reversals  of  Cyclic- 
Hardening. 


(b)  Branches  AB,  BC  and  CD 
are  superimposed  on  their 
initial  points  A,  B.  and  C, 
respectively. 


(c)  Branches  AB  and  BC  are- 
translated  along  the  elastic- 
slope  to  match  nonlinear 
portions. 


Figure  1.  Comparisons  of  hysteresis  branches  under  cyclic  hardening  phenomenon. 


linear  elastic  portion  is  termed  “yield  range  incre- 
ment" (YR1),  since  it  represents  a change  in  the 
yield  strength.  The  invariant  nonlinear  inelastic- 
portion  represents  an  intrinsic  material  stress-strain 
characteristic  that  is  independent  of  cyclic  history. 

Saturation  Behavior 

Anv  two  consecutive  hysteresis  branches  form  a 
hysteresis  loop.  Under  transient  conditions,  a 
hysteresis  loop  is  not  fully  closed  (Figure  la)  since  its 
two  branches  are  not  identical.  However,  in  the 
saturated  state,  a hysteresis  loop  is  fully  closed  since 
the  two  hysteresis  branches  are  identical,  which 
means  the  yield  strengths  or  YRls  in  both  branches 
are  equal.  In  most  metals,  the  saturation  value  of 
YRI  is  a function  of  the  stress  or  strain  amplitude. 
Based  on  available  data,  and  for  the  purpose  of  engi- 
neering analysis,  it  is  reasonable  to  assume  that  the 
saturation  YRI  is  independent  of  prior  cyclic  history. 
Microstructural  studies26  have  shown  that  certain 
types  of  materials  develop  a substruetural  disloca- 
tion cell  si/e  at  saturation  that  is  independent  of 
prior  history.  This  fact  lends  some  support  to  the 
assumption  of  prior  history  independence  of  satura- 
tion YRI. 


ni  1 belter  ami  ( I airil  "Cyclic  Stress  Strain  Response  of 
I t t Metals  anti  Alloys.-'  ACIA  Mi’tullurnu'a.  Vol  15  (October 
l%7). 


Figure  2 shows  the  saturation  YRI  plotted  as  a 
function  of  strain  amplitude  for  three  materials 
exhibiting  different  cyclic  characteristics.  The 
saturation  YRI  could  also  be  expressed  as  a function 
of  stress  amplitude  or  plastic  strain  amplitude.  In 
the  case  of  SAF  1045  normalized  steel,  which 
exhibits  cyclic  softening  at  low  strains  (approxi- 
mately <0.004)  and  cyclic  hardening  at  high  strains, 
the  saturation  YRI  increases  with  strain  amplitude. 
In  the  case  of  SAF  1045  quenched  and  tempered 
steel,  which  cyclically  softens,  the  saturation  YRI 
decreases  with  strain  amplitude.  Cyclically  harden- 
ing 2024-T4  aluminum  is  one  of  the  few  materials 
whose  saturation  YRI  is  independent  of  strain 
amplitude.  Figure  2 shows  some  evidence  of  prior 
history  independence  of  saturation  YRI  in  the  cases 
of  SAF.  1045  normalized  steel  and  2024-T4 
aluminum. 

Cyclic  Uniaxial  Properties 

Three  uniaxial  material  properties  appropriate  to 
the  cyclic  saturation  state  are  derived  from  a set  of 
hysteresis  stress-strain  loops  corresponding  to  that 
state:  cyclic  stress-strain  curve;  skeleton  stress-strain 
curve;  and  saturation  yield  strength  increment.  Such 
loops  may  be  obtained  from  fully  reversed  constant 
strain  cycling  tests  on  multiple  identical  specimens, 
each  specimen  being  cycled  at  a specified  strain 
amplitude  until  saturation  is  achieved.  Since  an 
ideal  state  of  saturation  may  never  be  achieved  in  an 


10 


SATURATION  rRi.KS'  M^o 


SOC554 


□ 

□ D □ 


O 

O 


Cvelie  Stress-Strain  Cunv 


7 


40C  '51 


n 

o 


□ 

o 


35(24  Ml 
3012071 

20(138) 


10(69)1 


4 14  ALUMINUM 


O CONSTANT  STRAIN 
HISTORY 

□ MIXED  STRAIN 
HISTORY 


SAE  1045  STEEL(NORMALIZEO) 

20(138) 

9 YR|»0. CORRESPONDS  TO 

\ MINIMUM  YIELD  STRENGTH 

\ CONDITION 

10(69) 

\SAE  1045  STEELlOfT  720*F) 

0 0 004  0 008  0 012  0 016  0 02 


STRAIN  AMPLITUDE 

Figure  2.  Saturation  yield  range  increment  for  three 
materials. 

actual  test,  the  condition  at  which  the  transient 
changes  are  minimal  may  be  considered  to  represent 
saturation  for  the  present  purpose.  Alternatively,  the 
saturated  hysteresis  loops  may  also  be  obtained  in  a 
single  multi-block  deeremental  step  test  as  illus- 
trated in  Figure  .1.  In  this  test,  the  specimen  is 
strain-cycled  at  a constant  amplitude  until  satura- 
tion is  achieved  at  that  level:  the  sample  is  then 
tested  at  the  nest  level,  and  so  on. 


This  curve  is  defined  as  the  locus  of  the  tips  ol 
lulls  reversed  saturated  or  stable  hysteresis  loops  ol 
various  sizes  (see  Figure  4).  Flic  hysteresis  loops 
shown  in  the  figure  arc  for  A -3b  steel  obtained  from 
a multiblock  deeremental  step  test.  An  altcrnatnc 
rapid  method  of  determining  the  cyclic  stress-strain 
curve  bv  means  of  an  “incremental  step"  test  is 
described  by  l.andgraf.  et  al.27 

If  the  hysteresis  loops  are  superimposed  on  their 
lower  tips  (Figure  5).  the  locus  of  the  upper  tips 
describes  a curve  that  is  geometrically  similar  to  the 
cyclic  stress-strain  curve,  but  magnified  by  a scale 
factor  of  two.  This  relationship  results  from  the 
symmetry  property  of  the  hysteresis  branches. 

Skeleton  Stress-Strain  Curve 

The  skeleton  stress-strain  curve  is  that  part  of  the 
stress-strain  characteristic  that  is  history-independ- 
ent; it  is  therefore  geometrically  similar  to  the  invari- 
ant nonlinear  portion  of  the  hysteresis  curve,  but  re- 
duced by  a factor  of  two.  A convenient  method  of 
determining  the  history-independent  portion  of  the 
hysteresis  curve  is  to  translate  the  hysteresis  loops 
plotted  in  Figure  5 along  the  elastic  slope  until  their 
upper  branches  are  matched  as  well  as  possible 
(Figure  b).  The  portion  OAB  of  the  hysteresis  curve 
in  Figure  b.  when  reduced  by  a factor  of  two.  will 
define  the  skeleton  stress-strain  curve.  I he  cyclic 
stress-strain  curve  and  the  skeleton  stress-strain 
curve  can  also  be  expressed  as  stress-plastic  strain 
relationships  by  eliminating  elastic  strain  compo- 
nents from  total  strains. 

Sn  In  rat  ion  Yield  Strength  Ineremenl 

The  saturation  Yield  Strength  Ineremenl  (YSI). 
which  is  one-half  of  YKI.  is  usually  expressed  as  a 
function  of  stress,  strain,  or  plastic  strain  amplitude, 
b is  given  by  the  difference  in  stress  level  between  the 
cyclic  stress-strain  curve  (Figure  7)  and  the  skeleton 
stress-strain  curve  expressed  as  a function  of  plastic 
strain.  As  the  name  denotes,  it  represents  the  change 
in  the  material's  yield  strength  as  a function  of  loop 
size. 


,TR.  W liimlgnit.  J Morrow. and  I l-.ndo.  ''Dclrrmination  ol 
I he  Cyclic  Stress  Strain  Curse.-'  ASTM  Journal  ol  Materials 
JMI  SA.  Vol  4.  No  I (American  Societs  for  testing  ami  Matcn 
als.  March  1%^).  pp  l7h  188 


Figure.!.  Strain-time  history  in  a multi-block 
deeremental  step  lest. 
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li  the  saturation  YSI  or  YR1  is  independent  of  the 
hysteresis  loop  si/e.  as  in  the  ease  of  2024-T4  alumi- 
ntint,  then  the  cyclic  stress  strain  curve  and  the 
skeleton  stress-strain  curve  differ  only  in  their  elastic 
portions  In  a constant  value.  These  materials  are 
designated  as  "Masing"  type  materials,  since 
Masing's  postulation  implies  such  a property.28 
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MULTIAXIAL  THEORY 
OF  CYCLIC  PLASTICITY 


Generalization  of  Uniaxial  Cyclic  Properties 


in  which  n = 0,  1 , 2 N - I . where  a von  Mises 

yield  condition  (given  by  f(0>  =0  for  plastic  tlow)  is 
used,  and 

p|j"  = o„  a|j"  \ d„(okk  u'k"k' ) |1  q 2| 

where  o >'*  = current  radius  of  the  surface  I’" 
oij  = current  stress  state 
djj  = Kronecker  delta. 

Fields  of  Plastic  Tangent  Moduli 
and  Isotropic  Hardening  Moduli 


T he  previously  discussed  observations  for  uniaxial 
cyclic  behavior  can  be  generalized  for  multiaxial 
stress  space  following  the  theory  proposed  by  Mroz29 
and  Iwan.30  The  skeleton  stress-plastic  strain  curve 
OS  shown  schematically  in  Figure  8(a)  is  approxi- 
mated by  N linear  segments  OA.  AB,  BC MN. 

Nodal  stress  values  corresponding  to  points  A,  B,  C. 

. . ..  N are  o‘°\  o‘u.  o12’ o'N  0,  respectively. 

In  a general  stress  space,  this  approximation  is 
represented  bv  nesting  hypersurfaces  f<0>,  fll),  f<2>, 
....  fiN)  (Figure  Mb)  with  initial  radii  o10>,  o'1’.  o<2). 

..  o (N  It.  respectively.  Initially,  these  hypersur- 
faces  arc  concentric,  with  their  centers  at  a point  0 
representing  zero  stress  level  in  the  stress  space.  The 
first  (smallest)  ot  these  surfaces,  f(0>.  is  the  conven- 
tional yield  surface  and  represents  elastic  limit. 

T he  hypersurlaces  translate  and  expand  or  con- 
tract during  plastic  flow.  Instantaneous  positions  of 
these  surfaces  in  the  stress  space  are  given  by  coordi- 
nates of  their  origins.  a|"'.  n = l),  I (N  — I). 

The  hypersurlaces  can  be  mathematically  repre- 
sented its 

f"”  = ;p;i'”pil;’l-'(o,,”)2=o  i Eq  1 1 
(summation  is  implied) 

)#H  K.  Jhansalc.  A Inctnm  Stress  Method  tor  the  C yclic  lit 
elastic  Behavior  ot  Metals,"  Iransai  Itont  ot  the  trit  International 
( ‘onlerenee  on  Strut  tural  Mn  harm  s in  Hear  tor  leehnolog  t Vol 
s l'.t|vr  1 s 4 i l‘rsi  . i ru  1 II  K Jhatts.ile  "A  Nc»  I’arametcr  lor 
the  Hssteretu  Stress  Strain  Hch.o  tor  ot  Metals,  Journal  of  l.ngi 
neermg  Material % ami  leehnologt  Vol  4"1.  M,,  | tlV’.S).  pp 
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A slope  at  any  point  of  the  skeleton  curve  of  Figure 
8a  is  defined  as  plastic  tangent  modulus.  This  defini- 
tion of  plastic  tangent  modulus  is  generalized  for  a 
multiaxial  stress  state  as  follows.  It  is  assumed  that 
the  projection  of  increment  of  stress  vector  dop  on 
the  exterior  normal  to  the  yield  surface  is  propor- 
tional to  the  plastic  strain  vector  dtfj.  T herefore,  if 
tiie  projection  of  doij  on  the  exterior  normal  to  the 
yield  surface  is  edt|’.  where  c is  a proportionality 
constant,  then  a vector  (do jj  - cdtj,)  is  perpendicular 
to  the  exterior  normal.  Mathematically,  this  can  be 
expressed  as 

,1(0) 

(doij  edej’)  =0  |Eq  T| 

'i  J o i j 

where  )f‘0,'Jojj  is  the  exterior  normal  to  the  yield 
surface,  and  the  proportionality  constant  c is  plastic 
tangent  modulus. 

In  Figure  8a.  the  skeleton  curve  of  the  uniaxial 
case  has  been  approximated  by  piecewise  linear  seg 

ments.  Nodal  points  0.  1 . 2 (N  I)  tire  at 

stress  levels  a101,  o'11,  o12’ olN  ",  respec- 

tively. I ,tch  nodal  point  m is  associated  with  a 
plastic  tangent  modulus  elm*.  which  remains  con 
stain  between  the  nodal  points  m and  tn  + I.  For  the 
multiaxial  stress  state,  this  approximation  is  genet 
ali/ed  bv  assuming  that  plastic  tangent  modulus  c ol 
I c|  .1  remains  constant  between  any  two  surfaces. 
Ibis  constant  value  of  c between  the  surfaces  fmit 
and  I ,m  * 11  is  given  bv  2c(nl1  4 where  e,ml  is  the 
corresponding  plastic  tangent  modulus  ol  the  piece 
wise  linear  skeleton  curve.  I Inis,  a field  of  plastic 
tangent  moduli  is  constituted  bv  associating  a plastic 
tangent  modulus  ctml  to  the  hypersurface  f<mi. 

\ field  ot  isotropic  hardening  moduli  is  defined 
similarly  as  the  field  ol  plastic  tangent  moduli  In  the 
uniaxial  ease,  stress  levels  o"  associated  with  the 
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PLASTIC  STRAIN 

(a)  Skeleton  stress-plastic  strain  curve. 


(b)  Hypersurfaecs  in  stress  space, 
t iKure  8.  Schematic  representation  of  hypersurfaces. 
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piecewise  approximation  of  the  skeleton  curve,  arc 
given  increments  of  do*'*  to  reflect  isotropic  harden- 
ing. All  these  increments  are  assumed  to  he  the  same 
to  insure  the  history-independent  nature  of  the 
skeleton  curve.  These  increments  are  assumed  to  he 
proportional  to  the  increment  of  plastic  strain. 

doth  = kdip.  i =0,  1.2 (N  - I)  |i:q  4| 

where  k is  the  constant  of  proportionality.  The  con- 
stant k ts  obtained  Iron i the  slope  of  the  saturation 
yield  strength  increment  versus  plastic  strain  ampli- 
tude curve. 

Tor  a multiaxial  stress  state.  Hq  4 is  generalized  as 

do1'*  = hdtp.  i = 0,  1.2 (N-l)  (Eq5| 

where  do*'1  = increments  in  the  radii  of  surfaces  f*'i 
h = isotropic  hardening  coefficient 
dt  P = equivalent  plastic  strain  defined  hy 


where  dr}’  is  increment  of  plastic  strain. 


Figure  V.  Piecewise  linear  approximation  of  yield 
strength  increment  versus  plastic  strain  amplitude 
curve. 

modulus  h,0>  is  chosen  to  approximately  reflect  the 
non  Maxing  and  transient  characteristics  of  real 
materials  in  the  uniaxial  ease: 


The  isotropic  hardening  coefficient  h is  assumed 
to  he  constant  between  any  two  hypersurfaces  f|m) 
and  f*m+  11  and  is  given  by  a value  h(n,)  associated 
with  the  surface  ftml.  This  constant  value  hlm)  is 
obtained  from  a piecewise  linear  approximation  of 
the  saturation  YS1  versus  plastic  strain  amplitude 

curve  (Figure  4).  Nodal  points  0,  I.  2 (N  1) 

on  this  curve  are  at  the  same  stress  levels  as  the  nodal 
points  of  the  piecewise  linear  skeleton  curve  and 

therefore  correspond  to  surfaces  f<0’.  f ‘ 1 ’ . f”’ 

ftN  li  respectively.  The  slope  of  a linear  segment 
between  nodal  points  i and  (i  + I)  is  designated  as 
k<*l. 

The  first  time  a stress  path  extends  into  the  zone 
between  surfaces  fh>  and  I'*'  * 1 > . the  isotropic  hard- 
ening modulus  hu)  corresponding  to  the  surface  f *'* is 
given  by  k *'* . For  any  subsequent  stress  excursion 
into  this  zone.  It1'1  is  set  equal  to  one-half  of  k*'* . ex- 
cept when  i 0.  This  factor  of  one-half  is  necessary 
to  properly  reflect  the  scale  factor  of  two  that  is 
associated  between  the  monotonic  and  cyclic  paths 
m the  uniaxial  case  (see  Cyclic  Uniaxial  Properties 
section  in  Chapter  4). 

f in  subsequent  stress  excursions  into  the  zone  be- 
tween surfaces  f‘01  and  I'*",  the  isotropic  hardening 


|do%  do”1 1 
h“»  = t y y 
AtT"1 


|Fq7| 


w here 


dov  = current  YSI 

do”’ and  At  I’*1’  = saturation  YSI  and  plastic  strain 
range  corresponding  to  node  I 
I = softening  hardening  factor  ~ O H 
to  1 . 

The  value  of  I I eliminates  transient  behavior  ti  e., 
a gradual  approach  toward  saturation);  I < I pro 
vides  lor  the  transient  behavior,  but  the  exact  satura 
lion  value  of  SSI  is  not  reached  A value  oil  OH  to 
0.0  seems  to  be  a good  compromise 

Translations  of  the  Hypersurfaces 

Translations  of  these  hypersurfaces  are  governed 
by  the  conditions  that  the  surfaces  do  not  intersect 
but  consecutively  contact  and  push  each  other  I hcv 
may  contact  only  at  a point  which  corresponds  to  the 
current  stress  state.  For  convex  and  geometrically 
similar  surfaces,  these  conditions  m.iv  be  satisfied  In 
requiring  that  instantaneous  motions  of  the  surfaces 
be  in  a specific  direction  This  direction  is  along  a 
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Figure  10.  Schematic  translations  of  hypersurfaces  in  stress  space. 


vector  connecting  the  point  of  contact  (current  stress 
state)  with  a point  on  the  next  surface  to  be  encoun- 
tered. which  has  a same  normal  direction  as  the 
normal  to  the  yield  surface  at  the  point  of  contact. 


— (i  + 1 ) 

0i+,9=  5,0,  <°,j  <’> 

Coordinates  of  point  0 are  given  by 


IN  01 


Translation  vectors  u ,,  . m =0.  I (N  — 1) 

can  now  be  obtained  by  considering  the  example 

given  in  Figure  10.  Surfaces  I'01.  tu> I "'are  in 

contact  at  a point  I’  corresponding  to  the  current 
stress  state  oir  Current  origins  of  these  surfaces  are 

denoted  In  O0.  0, 0,.  The  next  surface  to  be 

encountered  is  I"*  •>  with  origin  at  0i+).  Instan- 
taneous motions  ol  all  the  surfaces  will  be  along  the 
vector  1*0.  where  (,)  is  a point  having  the  same 
normal  to  f"  * 11  as  the  normal  to  the  yield  surface  at 
I’.  Stress  state  at  Q is  designated  as  projected  stress 
state.  Its  coordinates  arc  denoted  by  oH. 

To  obtain  coordinates  ojj.  a vector  0,  4*  Q is  draw  n 
parallel  to  ()„!’.  Homogeneity  of  F.q  I yields 

- ii -*  ti 

0,  . |0=<’o,0)  IN  HI 

where  o'0’  - current  radius  of  the  yield  surface 
and 

o"  * 11  current  radius  of  the  surface  I"  * ", 
I herelorc. 


Substituting  for  0j+  jQ  from  F.q  d. 


oH=„'i  + " + 5_(+0),,(o  «•••> 

'I  ij  o<0’  ')  " 


Instantaneous  motion  of  the  yield  surface  van  now 
be  written  as 


da'.®’ =dp(o‘j-ojj)  |Eq  HI 

where  da10’  is  an  increment  of  the  translation  vector 
ij 

ol  yield  surface  and  dp  is  a scalar  parameter.  F.q  1 1 
was  first  introduced  In  Mro/31  and  is  often  relerred 
to  as  Mro/'s  kinematic  hardening  law. 


The  parameter  dp  is  obtained  In  xnlisfving  the 
condition  that  the  current  stress  stale  remains  on  the 
yield  surface.  Appendix  A provides  the  definition  ol 
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an  expression  tor  dp  in  terms  of  stresses  and  stress 
increments. 


Coordinates  of  the  new  origin  of  the  yield  surface 
are  obtained  by  adding  du',j"  to  the  coordinates  of 
previous  origin.  Coordinates  of  new  origins  of  the  re 

maining  surfaces  in  contact  at  l*(f 1 I<J) f 1 1 ') 

are  delermineil  by  requiring  that  these  surfaces  re- 
main in  contact  at  the  new  stress  state  and  that  their 
origins  remain  on  a straight  line  joining  the  new 
stress  state  and  new  origin  of  the  yield  surface. 
Hence,  it  can  be  easily  shown  that 


,(m) 


7*111) 


= O 


to,, 


a10'). 

'J 


IHq  12| 


in  w hich  m = 1 , 2, 


The  surfaces  which  are  not  in  contact  at  the  cur- 
rent stress  state  are  also  translated  along  PQ  to 
retain  the  shape  of  the  skeleton  curve  under  uniaxial 
cyclic  loading.  Any  increments  in  the  radii  of  these 
surfaces  tend  to  distort  the  skeleton  curve.  To  offset 
this  effect,  origins  of  these  surfaces  are  translated  so 
that  measures  of  their  translations  are  equal  and 
opposite  to  the  increments  in  their  radii.  Therefore. 


da!m)  = dvrj  ,m=i  + l,i  + 2 N— 1 [Eq  13] 

i.i  'u  ^ 


where  dy  is  the  measure  (magnitude)  of  dajl[?l: 

dy  = — h('ld£p,  | Hq  1 4 1 


and  rj , , is  a unit  vector  along  I’Q: 


(o,|  0,|) 


' H°mn  ®mn^°nin  °ninH 


Hqs  5,  II.  12.  and  13  completely  specify  expan- 
sions (or  contractions)  and  translations  of  these  sur- 
faces. flic  plasticity  formulation  can  now  be  com- 
pleted with  the  following  How  rule”: 


d^dAf” 
" 8o„ 


IHq  151 


where  dA  is  a scalar  parameter  given  by 


dA 


»Td0M 


(2c  3)lol0,)J 


IHq  lh| 


Mill.  I he  Mathematical  I henry  of  Elasticity  (Oxford. 


Appendix  B provides  the  derivation  ot  Hq  15  and 
the  other  equations  used  in  the  plastic  analysis. 


Algorithm  for  Numerical  Analysis 


Ibis  section  gives  an  algorithm  for  calculating  a 


new  stress  stale  to  ),,,,,  at  the  end  ol  a given  strain 


increment  dr,,  by  (he  plasticity  formulation  pre- 
sented in  this  report.  The  following  variables  are 
assumed  to  be  known  at  the  beginning  of  the  given 
strain  increment: 


1 . Current  strain  i 


2.  Current  stress  o,, 


3.  ISURF  = number  of  the  largest  surface  in  con- 
tact at  the  current  stress  state,  (if  the  current  stress 
state  is  within  the  yield  surface,  then  ISURF  is  set 
equal  to  —I.) 


It  is  assumed  that  strain  increments  dr,,  are  small 


enough  that  changes  in  the  plastic  tangent  modulus 
and  the  isotropic  hardening  modulus  are  small.  If 
dt,,  are  not  small,  then  they  are  subdivided  into 
smaller  increments  before  the  following  algorithm  is 
applied. 


Step  1.  Define  new  variables  t,,  and  de,,  by 


Tij  = V 


dc'.j=d£il 


Assuming  an  elastic  behavior,  calculate 
d t,,  from  strain  increment  de,,.  using  Fq 
B-7  of  Appendix  B. 


Step  2.  Calculate  a trial  stress  state  T : 


Tij  = Tii  + dT..r 


Step  3.  a.  If  the  yield  function  f<0)  = fl#)(T|, . a l,0) 


o<0’><  0,  then  the  assumption  of  elastic 
behavior  is  valid.  Therefore,  the  new 
stress  state  is  given  by 


,0u)m»  1 n - 


ISURF  is  updated  by  setting  it  equal  to 
I I his  completes  calculations  for  the 
current  strain  increment. 


b II  f'01  X)  and  ISURF  = I.  then  trail 
sitioti  from  elastic  to  plastic  behavior 
lakes  place,  (io  to  Step  b. 


IX 


c.  II  l10>  >0  ;iml  ISURF  ^ 0.  then  plastic 
loading  continues.  Go  to  Step  4. 

Step  4.  From  the  strain  increment  de,r  determine 
stress  increment  il ti(.  using  Fq  HI  I.  Ma- 
terial constant  c in  this  equation  is  given 
by  the  plastic  tangent  iihhIuIus  e,n).  where 
n = ISURF.  Calculate  a new  trial  state  T,, 
bv 


Ti.i  =Ti, +dTir 


and  go  to  Step  5. 

Step  5.  Determine  f,",(T;j.  a!!1’.  o(nl).  n = ISURF 

1 u 


a.  Ilf""  <0,  then  the  stress  state  remains 
within  the  surface  f(l,),  n = ISURF  + 1. 
Therefore,  the  trial  stress  state  is  the  new 
stress  state,  as  shown  bv 


^°ij\icw  'if 


Updating  the  radii  and  origins  of  the 
hypersurfaces  (see  Chapter  4)  will  com- 
plete calculations  for  the  current  strain  in- 


lure. first,  the  present  approach  uses  a skeleton 
curve  derived  from  uniaxial  cyclic  tests  instead  ol  a 
uniaxial  monotonic  stress  strain  curve  I lie  second 
major  difference  is  that  the  yield  strength  in  other 
formulations  is  assumed  to  be  a function  ol  a mono 
tonically  increasing  parameter.  Ibis  restricts  the 
material  to  harden  or  soften  monolnnic.dlv.  regard- 
less of  the  nature  of  cyclic  loading.  Therefore,  in 
contrast  to  the  present  formulation  (where  yield 
strength  increments  are  defined  from  a field  of 
isotropic  hardening  moduli),  other  formulations 
cannot  exhibit  transient  hardening  or  softening. 

The  plasticity  formulation  developed  in  this  report 
is  capable  of  closely  simulating  all  the  important 
cyclic  material  characteristics’  under  the  uniaxial 
stress  state.  However,  to  demonstrate  the  real  worth 
of  the  present  approach,  its  predictive  capabilities 
must  be  verified  under  multiaxial  stress  states. 
Although  a few  multiaxial  cyclic  stress-strain  data 
are  available  in  the  literature,  attendant  cyclic  uni- 
axial properties  for  the  same  material  have  not  been 
generated  or  reported.  Since  the  present  formulation 
requires  cyclic  uniaxial  material  properties  as  a base, 
its  predictive  capabilities  cannot  be  tested  for  such 
data.  A limited  experimental  program  to  test  the 
present  formulation  under  a biaxial  stress  state  is  in 
progress. 


b.  If  f""  >0.  go  to  Step  b. 

Step  h.  Determine  a variable  RA  TIO  such  that 


|""(t  + RATIO*  dr  . u"".  oln>)  =0. 
il  y ii 

where  n ISURF + 1. 


Redefine  r1(  as  the  stress  on  f*n\  n = 


ISURF  4 I. 


T„  = T„+  RA  TIO  *dTj| 


Redefine  strain  increment: 


However,  significant  differences  between  the 
present  approach  and  alternative  approaches  cur- 
rently in  use  for  general  structural  analysis  (isotropic 
hardening  and  kinematic  hardening)  can  be  effec- 
tively demonstrated  lor  a typical  uniaxial  cyclic 
material  behavior.  The  simplest  and  most  important 
feature  in  a uniaxial  cyclic  response  is  the  saturated 
hysteresis  loop  under  constant  strain  limit  cycles, 
because  the  material  rapidly  approaches  this  state 
and  remains  there  for  a major  portion  ol  its  cyclic 
(fatigue)  lile.  Therefore,  the  accuracy  with  which 
this  saturation  hvslcrcsis  loop  is  predicted  gives  a 
reasonable  measure  ol  the  adequaev  ol  a plasticitv 
formulation. 


de,j  = ( I RATIO)  *deij 


Update  radii  and  origins  of  the  surfaces 
(sec  Chapter  4).  Update  ISURF  by  in- 
creasing it  by  one  and  go  back  to  Step  4. 

5 DISCUSSION 

There  are  two  major  differences  between  the 
present  formulation  and  most  others  in  the  litera - 


Figures  II.  12.  Id,  and  14  show  hysteresis  loops 
predicted  by  different  methods  for  two  materials 
A-dh  steel  (hot  worked)  and  SAT  1005  KtO'l  steel 
(cold  rolled).  Real  hysteresis  loops  ol  these  steels  (not 
shown  in  the  figures)  are  nearly  identical  to  those 
predicted  by  the  present  plasticitv  formulation  It 
the  plastic  tangent  modulus  ol  a material  is  non/cro 
lor  large  values  of  plastic  strain,  then  isotropic 
hardening  will  predict  continuouslv  growing  (or  con 
trading,  il  the  modulus  is  negative)  hvslcrcsis  loops 


1 


Figure  II.  Bilinear  hysteresis  loops  predicted  by 
isotropic  hardening  theory  lor  A-3b  steel. 

Figure  I I shows  such  hysteresis  loops  for  A-36  steel. 
It  the  plastic  tangent  modulus  approaches  zero  for 
large  plastic  strains,  then  hysteresis  loops  predicted 
by  the  isotropic  hardening  theory  will  stabilize  after 
a sufficient  number  of  cycles.  This  is  illustrated  in 
Figure  12  for  SAH  HXIS-HXW  steel.  However,  once 
these  loops  are  stabilized,  the  isotropic  hardening 
theory  reduces  to  a perfect  plasticity  theory  which 
does  not  allow  any  change  in  the  yield  range  of  the 
material.  In  the  case  of  the  kinematic  hardening 
theory  , yield  ranges  of  all  stabilized  hysteresis  loops 
remain  the  same  (twice  the  yield  strength)  and  slopes 
of  their  elasto  plastic  branches  are  equal  to  the  slope 
of  the  uniaxial  stress-strain  curve  at  a large  value  of 
plastic  strain.  Figures  Id  and  14  compare  the  stabi- 
lized hysteresis  loops  predicted  by  the  kinematic 
hardening  theory  and  the  present  plasticity  formula- 
tion. I hese  figures  clearly  show  that  stress-strain  re- 
sponses and  energy  absorption  estimates  (deter- 
mined Irom  areas  enclosed  by  the  hysteresis  loops) 
can  be  signiticantlv  erroneous  if  not  represented 
accurately.  I bis  error  becomes  more  critical  if  the 
material  is  subjected  to  irregular  or  variable  load 
histories. 
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APPLICATION  TO  SEISMIC 
DESIGN  PROCEDURES 


seismic  design  I oi  ev.iuipli  conxidci  the  tvpical 
scisntk  design  procedure  outlined  l>\  New  in. ok  and 
Hall"  and  shown  in  I iguie  I''  I lie  procedure  lias 
Inc  steps  t list,  an  c.uth>  kc  Ita/.iltls  oi  seismic 
nsk  is  selected  lot  a sue  based  on  the  geographic 
location  ol  the  slriiclure.  its  sicimtv  trout  active 
faults,  soil  conditions  at  the  site,  and  past  earth- 
quake activity  hi  that  region  Second,  a safety  lactor. 
allowable  deformation  limn,  and  allowable  prob- 
ability of  damage  are  selected  I hese  selections  are 
influenced  by  the  functional  requirements  of  the 
structure  and  the  costs  ot  a lailurc  in  terms  ot  human 
lives,  disrupted  services,  labor,  and  materials.  The 
values  selected  in  this  step  may  also  depend  on  the 
type  and  layout  of  the  structure  as  determined  in 
step  3.  Preliminary  design,  step  3.  is  accomplished  in 
accordance  with  design  codes  and  past  experience. 
Approximate  estimates  ot  the  structure’s  static  and 
dynamic  strengths,  ductility  , and  natural  frequency, 
damping,  and  energy  absorption  characteristics  are 
also  made.  In  step  4 the  overall  adequacy  of  the 
structure  is  verified  and  the  strength  or  other  param- 
eters are  changed  as  necessary . Steps  3 and  4 are  re- 
peated until  a satisfactory  design  of  the  structure  is 
achieved.  At  this  point  the  design  is  considered  final 
for  many  types  of  structures  for  w hich  a large  safety 
factor  is  implicit  in  the  design  criteria.  In  such  eases 
step  5 is  not  undertaken.  However,  it  is  important 
that  the  design  of  critical  facilities  such  as  hospitals, 
schools  and  nuclear  containment  vessels  be  verified 
by  an  analysis  (step  5)  in  which  structural  character- 
istics (strength,  damping,  energy  dissipation,  etc.) 
are  estimated  more  accurately  by  advanced  struc- 
tural techniques.  Such  an  analysis  may  indicate  need 
for  further  changes  in  the  design  and  a requirement 
to  repeat  steps  3.  4,  and  5. 

Energy  dissipation  characteristics  of  a structure 
are  estimated  (approximately  in  step  3 and  more 
accurately  in  step  5)  to  ensure  that  the  structure  can 
absorb  the  energy  imparted  to  it  In  ground  motion, 
blasi  forces,  or  strong  winds.  An  inahilitv  to  absorb 
this  energy  without  excessive  deformations  or  a long 
transient  response  may  result  in  catastrophic  failure. 
Energy  dissipation  is  usually  roughly  estimated  In 
assigning  some  damping  coefficients  to  the  si  rue 
lure.  These  coefficients  are  chosen  on  the  basis  ot 
past  experience  and  can  be  verified  onlv  bv  accur- 
ately determining  the  deformation  response  and 
associated  energy  dissipation  ot  several  tvpical  strm 
lures  Therefore,  the  research  conducted  under  this 


The  plasticity  formulation  described  in  this  report 
can  be  used  advantageously  in  several  areas  ot 


M Sew  mark  .iml  V\  J Mali  Pioirduivs  amt  t ntrru  lor 
I .tit lujii.ikt  KimnI.hu  IVsitfit  HuiUiini:  l*rtulut\  for 
^lut^nlton  MSS  iVitii»n.fl  Him. hi  .<1  St.nnl.ntU 
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Figure  12.  Hysteresis  loops  ot  SAL  1005  I OOd  steel  leolil  rolled). 


Figure  14.  Stabilized  hysteresis  loops  of  A-Ah  steel  (hot  worked). 


projeet  has  an  impact  on  seismic  design  procedures 
in  the  following  areas: 

I.  Results  of  accurate  inelastic  analyses  (step  5) 
will  preside  guidelines  for  developing  structural 
strength  criteria  tslcp  .1)  under  inelastic  loading  con- 
ditions. 

I nergs  dissipation  caused  In  cyclic  inelastic 
deformations  will  provide  a rational  basis  for  deter- 
mining damping  coefficients  at  high  stress  levels  (to 
be  used  in  step  A). 


.V  Conservative  and  elastic  designs  will  be  re 
placed  by  inelastic  designs,  resulting  in  substantial 
savings  in  labor  and  materials  costs. 

4.  Cyclic  inelastic  analyses  will  accurately  deter 
mine  the  deformation  response  and  safety  ol  existing 
critical  structures  (step  5). 

5.  Ihe  capability  to  perform  an  accurate  strttc 
tural  analysis  will  enable  designers  to  use  novel 
design  concepts  with  confidence.  Ihis  will  allow 
more  creativity  in  the  use  of  space  and  materials 


’.1 


Figure  15.  A seismic  design  procedure. 


7 SUMMARY  AND  FUTURE  WORK 
Summary 

This  report  has  presented  a plasticity  formulation 
adequate  for  cyclic  inelastic  structural  analysis; 
appropriate  cyclic  material  properties  required  by 
this  formulation  can  be  determined  from  simple 
cyclic  uniaxial  tests.  This  formulation  is  capable  of 
simulating  all  the  important  cyclic  material 
phenomena  of  deformation  response,  including 
memory  of  prior  history  and  non-Masing  type  of 
behavior.  The  need  for  classifying  a general  loading 
path  into  a plastic  loading  or  a plastic  load  reversal 
cyclic  analysis  has  been  eliminated.  Hysteresis  loops 
predicted  by  the  present  formulation  are  very  close  to 
a material's  actual  hysteresis  loops.  These  loops  have 


been  shown  to  be  significantly  different  from  those 
predicted  by  other  plasticity  theories  currently  used 
in  finite  element  codes.  Therefore,  using  the  present 
formulation  in  existing  nonlinear  programs  will 
achieve  a reliable  and  accurate  estimation  of  energy 
dissipation  and  deformation  response. 

Future  Work 

Future  research  tor  this  project  includes  tl)  ex- 
perimental verification  of  the  plasticity  formulation 
presented  in  this  report  under  a cyclic  tension- 
torsion  state,  and  (2)  an  analysis  of  multiaxial  in- 
elastic hvsterctic  behavior  of  reinforced  concrete 
structural  members.  The  results  obtained  in  this 
project  will  be  translated  into  simplified  design  pro 
ced ures  to  be  incorporated  into  the  Seismic  Design 
nl  HiiilJinfi\.  I Ms  HI  to  |()  series. 
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APPENDIX  A: 


DERIVATION  OF  AN  EXPRESSION 
FOR  THE  PARAMETER  dw 


Partial  differentiations  of  fl0>  with  respect 
u‘0’  are  given  by 


The  current  stress  state  remains  on  the  yield  sur- 
face during  plastic  flow;  therefore, 

df<0>  =0  [Eq  All 


where  f<0)  is  given  by  Eq  I . 


3f>»>  = _3f(0> 

do, | da'*’  >J 

From  Eqs  A2,  A3,  and  1 1 

P^>  dojj  - P‘«>  dM(oH  -o|()  3<(”  d5'#’ 


Eq  A I yields 

H'to>  jfio)  ,r<0) 

^do  +,*c.)da»»+  ,^ydo<o>=0  |Eq  A2| 


Therefore. 


C,  domn  +j_5‘°’d-o‘«- 
Pl0>  (o?.-o  ) 

U IJ  tj 


to  o,,  and 
I Eq  A3| 

= 0. 

I Eq  A*»| 
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APPENDIX  B: 


DERIVATION  OF  EQUATIONS 
OF  INCREMENTAL  PLASTICITY 

From  Eq  1,  yield  surface  is  given  by 

f(o»  = p.0,  p«o,  _ 2 (5<0))2=o  |EqB1| 

Differentiation  of  Eq  B1  with  respect  to  o gives 


= put, 
30j|  'i 


|Eq  B2| 


From  Eqs  15  and  B 2.  the  flow  rule  can  be  written  as 
dr!;=dAP;o>  |EqB3] 

Substitution  of  Eqs  B2  and  B3  into  Eq  3 yields 


(dojj  — cdA  P‘?>)  P‘°>  =0. 


Therefore. 


dA  = 


..  pio)  p(b) 
0 1 1 


From  Eqs  B1  and  B4 


I Eq  B4| 


dA  = 


Pfd°U 

(2c  3)(o<0’)2 


(Eq  B5| 


fo  obtain  dA  in  terms  of  strain  increments  dt,,. 
total  strain  increments  are  decomposed  into  their 
elastic  and  plastic  components 


d«u“d,e+d,B 


I Eq  BO  | 


From  elastic  constitutive  equations,  do 


as 


is  v\  ritten 


do.i  = Eijkt  d«ki  IFq  B7| 

where  E,,k|  are  elastic  coefficients. 

From  Eqs  B(>  and  B7. 


do.,  = E,,kl  <d«kt-d4V  (Hq  BK| 
Eqs  B3  and  BK  yield 

doM  = H,,kt  <d£k,-dA  Pff)  (Eq  B4| 

By  substituting  do,,  front  Eq  B4  into  Eq  B5  and  re- 
arranging terms.  dA  can  be  written  as 


dA 


fTE.,k,  d‘k, 

(2c  3)(o‘0>)  + P<0)E  Pl0> 
r\  rstu  1 in 


I Eq  BIO) 


dA  from  Eq  BIO  can  now  be  substituted  into  Eq  B4 
to  express  do,,  in  terms  of  total  strain  components 
d£kf 


do,,  = 


JXO)  L r pi 0) 

nil)  “mnkl  i||>q  pq 

■qkl  (2c  ^)(o(0>)2  +|>(0)|  |*<0> 

rs  'rslir  in 


I di  - 


II  q Bill 

Material  constant  c used  in  Eqs  BIO  and  Bll  ts 
obtained  from  a field  of  plastic  tangent  moduli  as 
explained  in  Chapter  4. 
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